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ABSTRACT 


We consider three-neutrino MSW mixing, assuming m3 m2 > mj as expected from 
theoretical consideration if neutrinos have mass, and calculate the corresponding mixing 
parameter space allowed by both the 87 Cl and Kamiokande II experiments. ^^e also 
calculate the expected depletion for the 71 Ga experiment. Finally, we explore a range 
of theoretical uncertainty due to possible astrophysical effects by varying the 8 B neutrino 
flux and redoing the MSW mixing calculation. 
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I. Introduction: 

The well known solar neutrino problem centers on the fact that the measured time 
averaged solax neutrino flux is lower than theoretical estimates by a factor of 2~3 in two 
experiments, the Homestake 3 'Cl experiment and the Kamiokande neutrino-electron scat- 
tering experiment. In the 37 C1 experiment, the observed neutrino flux=2.2±0.3 SNU (Solax 
Neutrino Unit, 1 SNU=10 _36 capture events/target nucleus/sec) 1 , compaxed with the pre- 
dicted 7.9 SNU by the Standard Solar Model (SSM) of Bahcall et aJ 1 or 5.8 SNU of the 
French solar model. 2 Although many attempts have been made to assign an estimate of the 
uncertainty 3,4 in the theoretical numbers, the spread between the French and US models is 
probably as good as any estimate (it should be noted that the two models do approximately 
agree when the same choice of input parameters are made but differ primarily because of 
different choice as to the preferred input parameters). In the Kamiokande II experiment, 
the measured fluxes=0.45 ± 0.06 ata< ± 0.06 aj , at relative to the predicted value of Bahcall et 
a!. 5 The Kamiokande and Homestake experiments are mainly sensitive to the higher energy 
8 B neutrinos with the latter also have some admixture of 7 Be neutrinos. 1 Furthermore, 
two new experiments using 71 Ga (GALLEX and SAGE) which can see the main-line pp 
neutrinos have recently begun operating, while preliminary results from SAGE also appear 
reduced relative to theory. 6 These reports remain to be established following subsequent 
calibration studies using 51 Cr sources to verify that ' 1 Ge produced by neutrino capture on 
71 Ga can indeed be quantitively detected. 

Vacuum neutrino mixing was proposed to solve the problem, but it required large 
vacuum mixing angles relative to the Cabbibo angle ( 6 C = 13°) to provide a maximal 
flux reduction of 1 fN v , where N„ is number of neutrino flavors involved in mixing. 1 The 
Mikheyev-Smirnov-Wolfenstein (MSW) matter mixing effect 7 ’ 8 enables neutrinos to flip 
flavor while propogating in matter even with a small vacuum mixing angle. Thus, it pro- 
vides a more natural solution to the solar neutrino problem. 9 In the two-flavor mixing case, 
it has been found that in order to get the observed reduction with the MSW effect, only 
a triangular area in parameter space A = (m 2 - m\) vs.sin 2 20 12 fcos29 1 2 is allowed. 1 Full 
three-flavor MSW effects have also been investigated previously by Kuo and Pantaleone 10 
and Barger et ad 11 . In this paper we reexamine the full three-flavor mixing in the light of 
the present experimental situation. Since we now know from LEP 12 and from cosmology 13 
that there are indeed 3 neutrino families, it seems only reasonable that the formalism used 
should take into account this realty. However, we will see that the basic two-neutrino 
picture still gives a reasonable understanding of most of the action. We will also assume 
m 3 ^ m 2 > mi, which, if neutrinos have mass, is a reasonable theoretical expectation. In 
section II, we derive the probability of solar neutrinos flipping their flavors on the way to 
the earth. In section III we calculate numerically the allowed parameter space from the 
data of the 3 ‘ Cl and Kamiokande II experiments assuming three- flavor matter mixing. 
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Mixing between three othogonal neutrino states can be generally described by a unitary 
operator U 

K) =U Qi \ui), a = e,n,r, *'=1,2,3 (1) 

where u Q denotes flavor eigenstates and i y , mass eigenstates. 

Before looking at the complication of three-flavor mixing, let us review the standard 
two-flavor formalism. 14 ’ 15,16 In the special case of two neutrino mixing between v t and v x 
(v x might be a linear combination of and u T ), 

tt _ f cos6 sind\ 

U ~ \-sin6 cosd) ■ K ) 


The Schrodinger Eq. of a neutrino travelling through 


matter in flavor basis is 


. d 


i — 

dt 



( 3 ) 


Let A=ml — mf, 


_ _j_ f —Acos28 + 2A Asin2d\ _ f —cos28m sin28\f\ 
~ y Asin2d Acos26 J ~ 4 E \ sin2$M cos26m ) 


( 4 ) 


where A = 2 y/2G pn e E is the induced mass of the electron neutrino due to the MSW effect. 
Gp is the Fermi constant, n e is electron density of the matter and E the neutrino energy. 
A term proportional to the identity matrix has been subtracted. In the second expression, 
Am and are counterparts in matter of A and 8, where 


A m = \{Acos28 - A) 2 + (A sin28) 2 }^ 2 (5) 

tan28M = tan28/[l — (A/ A)stc28]. (6) 

When .4 = Acos29 , 9m = 7 r/ 4 , A m reaches its minimum, the two mass eigenstates are 
almost degenerated and the neutrino encounters a resonance. 

An interesting case occurs if neutrinos propagate through matter which has a slow 
and monotonically varying electron density. In that case the neutrinos can pass through 
the resonance adiabatically and flip their flavor (Fig.l). But generally, neutrinos evolve 
adiabatically in regions away from resonance. At resonance, they have some probability 
to jump between two mass eigenstates and thus don’t change their flavor completely. 
Such jump probability can be obtained by solving the Schrodinger eq.(3) assuming linear 
variation of electron density at the resonance region. The result is 1 7 


— 7rAsm 2 20 
Pjump = exp[ Tgr— 


n fi 


4Ecos26 y \dn e /dr\ 


) /*e3 < • 


( 7 ) 
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The probability for a i/ e produced at a high density region, passing through resonance, 
reaching the vacuum region, to maintain its flavor is 14 


P — 2 + ( 2 -Pj ump)*-O'S2$^/COs20 


(S) 


where 9 m is the mixing angle at which the u e is produced. If the density at which v t is 
produced is sufficiently large, 9 m —■ * tt/2 (see eq.(6)) and eq.(8) becomes 


P & sin 2 9 + Pj ump cos26. 


(S') 


To obtain the observed neutrino flux reduction, a large fraction of the i/ e ’s produced 
at the center of the Sun have to change their flavor. This requires the following: (a) 
part of the solar neutrinos have to pass through a resonance; (b) from eq.(8) and (8'), 
Pjump must be sufficiently small; (c) also from eq.(8) and (8') sin 2 # has to be sufficiently 
small. These conditions yield constraints on neutrino masses and mixing angles and limit 
them to a triangle in parameter space A vs.sin 7 29/cos29. 1 Each region of the triangle 
produces a different spectrum for the resultant neutrino flux to be observed. On the 
horizontal region, only high energy neutrinos go through resonace and change flavor; thus, 
the high energy part of neutrino spectrum gets “cut off”. This region is effectively ruled out 
because the Kamiokande II experiment which is only sensitive to high energy neutrinos 
(Ethrcakold of recoil electron=7.5MeV 5 ) detected a larger fraction of its theoretical flux 
than the 37 Cl experiment (Etft rej ,/ lo /d=0.814MeV 1 ) did. On the diagonal region, the lower 
energy neutrinos go through resonance with a smaller Pj ump , they are thereby supressed 
more than the higher energy ones. This is consistant with experiments. On the vertical part 
of the triangle the entire spectrum is uniformly suppressed. This is also not contradictary 
to the observation at the 3a level. 


II. Matter mixing for three-neutrino mixing: 

The three-neutrino matter mixing is more complicated. It involves three mixing angles. 
The unitary operator 

U = e^ A7 e -«^5 e iu,A 2 ( 9 ) 

where A 2 , A5, A7 are the generators of SU(3), 

/0 -i 0\ /0 0 -i\ (0 0 0 \ 

A 2 = i 0 0 ,A 5 = 0 0 0 ,A 7 = 0 0 -i . (10) 

\0 0 0 / \i 0 0 J \0 i 0 / 

In our problem, we are only interested in how much the v e flux is reduced without caring 
about which flavor it is converted to. We can rotate the flavor basis by exp(— to get 
rid of ip. The new basis is i/ e , 1 / p i, v r i and 

( cosocosu cos<p sin u —sin<f>\ 

—sinuj cosu 0 j . (11) 

sin<f>cosu sin<j)sinu> cos<f> J 
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In general there are two resonances. These resonances can be coupled so it is hard to see 
when they happen and how neutrinos behave while crossing these resonances. However, if 
7 Ti 3 » m 2 > mi , the resonance between two lower mass eigenstates can be well decoupled 
from the higher one. The assumption that m 3 m 2 > mi is reasonable from see- 
saw model consideration . 18 This also would be anticipated if v T , the heaviest neutrino, is 
associated with Hot Dark Matter where m„ r ~ 25eV. 19 Kuo and Pantaleone 10 found that 
under such assumption the lower resonance happens when 

Acos 2 <f> = Acos'2lv (12) 


and the higher resonance occurs when A = m\cos2<}>. We will follow the above assump- 
tion to solve the Schrodinger equation for three-neutrino mixing. We also assume that 
the higher resonance requires a much higher density than that of the Sun and becomes 
irrelavant. To find P jum p for the lower resonance, we further rotate the basis by exp(z<^A 5 ). 
The new basis becomes v e > , v T » (both linear combinations of u e and u T i) and The 
associated Schrodinger equation is 


4 ( < 4 $ | = M 



with 10 



—Acos2u> + 2 Acos 2 4> Asin2u> 
Asin2u Acos2u ? 

—Asin2(p 0 


—Asin 2<j> 

0 

2m\ — m\ — m] + 2 Asin 2 <j> , 


(13) 


(14) 


Under the assumption m 2 m 2 and m 2 , we can see that C T « oscillates so fast that 
it averages out (if we take the see-saw model, 7773/7712 ~ m 2 ./m 2 ll ~ 10 2 using leptons (or 
~ m 2 /m 2 ~ 10 4 using quarks), C T n would oscillate ~ 10 4 (10 8 ) times faster than the other 
two components). Then the transition probability of the lower resonance Pjump is given 
by eq.(7). 

If we let 

^1,2 or 3(^42) = exp [-7 ^ £1,2 or- 3(0*] > ( 15 ) 

E\,E 2 ,E$ are energy eigenvalues of three mass eigenstates, the time evolution of a v t 
produced at time t % , passing through lower resonance, detected at time tf is 


+ U e2 {ti)^ 2 (ti^t r )\u 2 (ti)} + t r ) 1 2^3 (^t )) 

(before resonance) (16) 


and 

- a2U e2 (ti)$2(U,tr)]Wl(t /)) 


5 


+ ‘$2(^r, t f)[a*U e 2(ti)$2(U, tr) + «2^el (t»)^l (*»'?*r)] I ^2 ( ^ / ) ) 

+$ 3 (/ r ,i/)a 3 ^ e 3 (<i)$ 3 (ii,< r )|y 3 (f/)) (after resonance) (16') 

where |a 2 p = Pj ump cos 2 <j) (remember the Pj U mp we calculated earlier was in the basis 
i/ e /,ty and u T n, a rotation of basis i v T > by angle <^) , |ai| 2 = 1 — Pj ump cos 2 <f>, and 
|a 3 p, from the above decoupling approximation, is 1. The amplitude of detecting a v t at 
tf is 

A e ,e = [aib T e j(ti)$i(ti,i/) — a,2U e 2(ti)$l(tr,tf)$2(ti,tr)]Uel(tf) 

+[a*U e2 (ti)$2(ii,tf) + a 2 U e i(ti)^ 1 (t i ,t r )^ 2 (t r ,tf)]U e 2(tf) 

-\-a 3 U e3 (ti')U c3 (t y)$ 3 (tj,ty). (17) 

After averaging over the position of production, i.e. <*, and the position of detection t/, 
the probability P of survival of v t is 

P = \aiU e i(ti)U e \(tf)\ 2 + \a* 2 U e2 (U)U el (t f )\ 2 + \a\U e2 {U)U t2 {t f )\ 2 

+\a 2 U el (t t )U e2 (t f )\ 2 + \a 3 U e3 (ti)U e3 (t f )\ 2 . ( 18 ) 

Taking expression of U ei from eq.(ll) 

C7 e i = cos4>cosu , U e 2 = cos<j>sinu, U e3 = —sin4> (19) 


we get 


P = [7 + (- — Pj UTnp cos 7 4>)cos2uJMCos2u>]cos 2 <l>MCOs 2 <f> + sin 2 4>M$in 2 4> 
2 2 


+ (i — Pj UTnp cos 2 4>)cos2u>MCos2u]cos 4 <j) + sin 4 <f> 


( 20 ) 


where <}>m and u>m are and u>(<,), the mixing angles at the position of production; 

4>{tf) and io(tf) are <p and u> if we detect in vacuum. We can see that the u e flux will suffer 
an additional overall supression when <f> ^ 0. 

Prom eq.(20) we can calculate the experimental constraints on parameter space A vs. 
sin 2 2uj / cos2u , which reduces to A vs. sin 2 26 / cos20 when <j> = 0. 


III. Numerical Result: 

Figures 2-4 show the allowed parameter space for different <f> with 1 <t and 3cr errors, 
from the results of both 3 'Cl and the Kamiokande II experiments (by convention, we 
replace u by 9). Fig. 2 (b) is basically consistent with Chen and Cherry’s result. 20 We can 
see that the horizontal region is still ruled out, in agreement with our previous discussion 
of the two-flavor mixing case. It is interesting to note that at the la level (Fig. 2(a) to 
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4(a)) there is no overlap for the vertical region either. This is understandable: for the 1 a 
error, the lower bound of the Kamiokande II is 0.45 — \/0.06 2 + 0.06 2 = 0.365, which is 
higher than the upper bound from 37 Cl experiment 2.5/T.9 0.32; besides, on the vertical 

region, the neutrino spectrum is uniformly suppressed. Thus, there can’t be any overlap. 
So at the la level the only possible parameter space is the diagonal region. 

From Fig. 2(b) to Fig.4(b) we find at the 3 a level that the vertical region is then 
allowed. However, part of this region will show day-night effect because of the regeneration 
of electron-neutrino flux in the earth. 21 So, after considering the Kamiokande IPs limit on 
the absence of a day-night effect, region sin 2 20/cos20 > 0.02 and A=2xl0 -6 -10 -5 eV 2 is 
also excluded at 90% C. L. when <j) = 0. 22 At large <f>, the overall suppression factors in 
eq.(20) which are not affected by neutrinos travelling through the earth (because of its 
low density) will smooth out the regeneration. Therefore we expect the excluded region 
to shrink. 

When <£=0.1rad (Fig. 3), we find no significant change in parameter space. This case 
would be interesting if the reported 1% mixing of 17keV neutrino in /? decay turn out to 
be true (then the 17keV neutrino would serve as the heaviest neutrino species). 23 At large 
4> , Fig. 4 shows that the allowed parameter space is shifted, but we still find the same 
situation with respect to the vertical and diagonal regions as in <j>=0 case. 

The ongoing 71 Ga experiments have a very low threshold 0.233MeV. Thus, they will 
bring more severe constraints on the parameter space. As mentioned earlier, at present 
SAGE gives a preliminary low result. 6 This would suggest that the central diagonal region 
would be the final answer if the MSW effect is the solution to the solar neutrino problem. 

In standard Solar Models, the 8 B neutrino flux is very sensitive to the temperature 
at the center of the Sun (a T c 18 ) whereas the other neutrino sources are less sensitive. 1 ’ 24 
Therefore the theoretical expected fluxes calculated for current 8 B dominated solar neu- 
trino experiments depend drastically on what T c is calculated in the Solar Model. This 
gives rise to greater astrophysical uncertainty for the determination of neutrino parameter 
space. To try to estimate the sensitivity to the uncertainty we have calculated the allowed 
parameter spaces for different central temperatures. If we assume the central temperature 
is only slightly varied, the predicted solar 8 B neutrino flux would be changed by some fac- 
tor while other neutrino fluxes are almost unaffected. Fig.5(a) and (b) show the overlap of 
two experimental results under different solar 8 B fluxes at la error. We see that a change 
of the solar 8 B flux by a full factor of 2 will only shift the overlaped region back and forth 
without inducing any new allowed parameter space (after the completion of this work, we 
notice that Smirnov also calculated the allowed parameter space for different T c , but at 2 a 
level 25 ). We also find that at la error, in order for the neutrino to have any allowed MSW 
parameter space, the predicted 8 B neutrino flux for the 37 Cl experiment must be higher 
than ~2.5 SNU when <j> = 0 and ~5 SNU when <j> = 0.5rad. For the 71 Ga experiments, 
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such changes in the the solar S B neutrino flux produces neglegible effect because Gallium 
experiments are dominated by the pp neutrino flux, v_ 

Based on Figures 2-5, we can conclude that if 71 Ga does give a very low flux, A = 
m 2 — m i w iU i n 10“ ' to 2 x 10 -6 eV 2 , and sin 2 20 will be 0.01-0.3, which is natur all y in 
agreement with Bahcall and Bethe’s prediction 26 (using only two neutrino families and a 
similar analysis of the 37 Cl and Kamiokande implications) and consistent with the Cabibbo 
angle 6 C ( sin 2 28 c = 0.2). This corresponds to a mass of ~ 10 -3 eV, which is more 
them 10 3 times larger than that of some SUSY see-saw models, and ~ 5 times smaller 
than the lowest m Ufi value in the SO(IO) see-saw model, both of which are discussed by 
Bludman et a l. 27 Thus our is within the admittedly large theoretically predicted range. 

Generally, if we take the see-saw inspired relationship = m\/m\ x C(C 0(1) 

but larger than l) 27 , using m c = 1.55GeV, m t = 124 ± 34GeV, 28 we obtain m Vr ~ lOeV. 

This would make v T a natural candidate for the Hot Dark Matter (see also discussion in 
ref. 29). Similarly, we obtain from the see-saw model m„ e ~ 10 -8 eV. These n um bers are 
consistent with the current experimental bounds on the neutrino masses and the neutrino 
oscillation parameters (Fig. 6). 30,31 

IV. Conclusion: 

After taking into account the three-neutrino matter mixing under the assumption 
that m 3 » m 2 > m ly we find that for various values of the second mixing angle <j>, 
the neutrino mixing parameters most probably lie in the diagonal region of parameter 
space Avs.sin 2 29 / cos26 . The vertical region is less probable but hasn’t been ruled out at 
3a level (except those parts ruled out by the lack of a day-night effect). The horizontal 
region of the parameter space is excluded by current experiments. The Gallium experiment 
may give us a final resolution between the diagonal area and vertical area. Changing the 
high energy 8 B neutrino flux will shift the allowed parameter space, but will not change 
the shape of allowed space significantly. The basic robustness of two-flavor MSW mixing 
results is not weakened by taking into account the full three neutrino flavors. 
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Figure Captions: 


Fig. 2-4 

The doted region is allowed by the 37 Cl experiment; the crossed region is allowed by the 
Kamiokande II. Percentages at the lower left show expected flux for the 71 Ga experiment 
relative to SSM. 

2(a) <£= 0, la error. 2(b) <£= 0, 3cr error. 

3(a) <£=0.1, la error. 3(b) <£=0.1, 3<r error. 

4(a) <£=0.5, la error. 4(b) <£=0.5, 3 a error. 

Fig.5 

This figure shows the parameter space allowed by both the 37 Cl experiment and 
Kamiokande II with different predicted solar 8 B neutrino fluxes for the 37 C1 experiment 
assumed. 

5(a) <£ = 0, la error. 5(b) <£ =0.5, 1<7 error. 

Fig.6 

The current experimental constraints on neutrino oscillation parameters. [1] — ♦ v e \ 

[2] Vp <-> i/ T ; [3] v t *-> v T . [4] u e — > u x ] [l], [2] and [3] are from accelerator experiments, [4] 
is from reactor experiments. The left side of the lines is excluded at 90% C. L.. 
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